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ABSTRACT 



We present results of nearly six years of spectroscopic observations of the B[e] star V2028 Cyg. The presence of the cold- 
type absorption lines combined with a hot-type spectrum indicate the binarity of this object. Since B[e] stars are em- 
bedded in an extended envelope, the usage of common stellar atmosphere models for the analysis is quite inappropriate. 
Therefore, we focus on the analysis of the long-term spectral line variations in order to determine the nature of this 
object. We present the time dependences of the equivalent width and radial velocities of the Ha line, [Ol] 6300 A, 
Fen 6427, 6433, and 6456 A lines. The bisector variations and line intensities are shown for the Ha line. The radial 
velocities are also measured for the absorption lines of the K component. No periodic variation is found. The observed 
data show correlations between the measured quantities, which can be used in future modelling. 
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1. Introduction 



binaries: general - circumstellar matter - Stars: emission-line, B[e] - Stars: individual: 



V2028 Cyg (MWC 623) belongs to a rather heterogeneous 
group called B[e] stars. This group includes such dissim- 
ilar objects as compact planetary nebulae, Herbig stars , 



of both the Balmer series and low-ionized metals originate. 
The polar areas are occupied by a CAK (Castor, Abbott 



and supergiants or symbiotic stars (Lamers et al. 1998) 
A large number of B[e] stars still belong to this "unciassi 
fied" subgroup, which is also the case for V2028 Cyg. It is 
difficult to assign this object to a specific s ubgroup because 
o f the presence of a circumstellar matter. [Miroshnic henko 
| ( |2007| found that many unclassified objects display com- 
mon spectral signatures. He denoted this group as FS CMa 
stars, since this star can be taken as a good prototype 
for these objects. He included in this group V2028 Cyg. 
Current observation analysis shows that they are non- 
supergiants. FS CMa objects are probably more evolved 
stars, but are either close to the main sequence, or stil l 
on it |Miroshnichenko et al. ||2000| | Miroshnichenko [2007 1. 
Miroshnichenko | ( |2007p suggested that a binarity could ex- 



& Klein 1975) wind of a low density and high velocity. 
The gas in this outflow is hot and the emission lines come 
from metals in higher ionization stages. 

V2028 Cyg was firs t noted as an emission- l ine star 



by |Merrill et aT~| ( |1942[ ). |Allen fc Swings | ( |1976| ) discov 



ered its photomctrical variability (amplitude 0.5 mag in fil 



ts p 

ter V). [ Bergner et aL 



( 1995[) the n presented the ampli- 
tude 0.2 — 0.4 mag. Allen ~ fl973 ) analysed photometric 
measurements from the filters H, K, and L and reported 
an infrare d exces s that he explained in terms of a dust 
envelope. 



Allen 



plain the observed properties of these stars. 

The commonly used m odel of B[e] stars was proposed 



(1974) then noted the connection be- 
tween the_excess and for bidde n emission lines. A few years 
later, Allen & Swings (19761 published a spectral analy- 
sis of the group of peculiar Be stars with infrared excesses, 
including V2028 Cyg. They identified permitted and for- 
bidden emission lines, such as Balmer lines, Fen, [Ol], 
and [Fe n] , which dominate the emission spectrum. 

V2028 Cyg was noted for t he first time as a binary 



by Zickgraf et al. ([1985) . The B[e] phenomenon arises 
in a two-component stellar wind with a dusty ring. A slow, 
cool, high density wind occurs around the equatorial plane. 
This wind forms an outflowing disc, where the emissions 



by Arkhipova fc Ipatov 
the types B8 and M1III 



Send offprint requests to: J. Polster, e-mail: jpolster@email.cz 
* Based on data from the Ondfejov 2 m telescope, Czech 
Republic. 



(11982). They suggested spectral 
or the system from the analy- 
sis of spectrophotometric observations. However, they used 
assumptions that were l ater found to be incorrect. High res- 
olution spectra enabled |Zickgraf fc Stahi~| ( |1989[ ) to iden- 
tify many absorption lines of mainly neutral met als, such as 
Fei, Til, Vi, Cai, an d Li I (also detec ted later by Corporon 
fc Lagrange |1999 and Zickgraf 2001 ). By fitting theoretical 



1 



Polster et al.: Time-dependent spectral-feature variations of V2028 Cyg 



spectral energy distribution to photometric data, they sug- 
gested that V2028 Cyg is a binary composed of a B2 dwarf 
and a K2 giant . From comparison with spectra of normal 



stars, 
to B4 



Zickgraf (2001) refined the component classification 
11 and K2Ib-II. Multico lour optical and near-i nfrared 
photometric observations led Bergner et al. (1995) to de- 
rive the type K7III for the cool component. The primary 
component falls on the HR diagra m into the region wher e 
other FS CMa stars are also found ( Miroshnichenko |2007 ) . 
Its position close to the region occupied by classical Be 
stars, provides a good opportunity to study the connection 
between the Be and B[e] phenomenon. 



Zickgraf (20011 derived a kinematical distance to- 
wards V2028 Cyg using the average radial velocity 
(2.0^Q 3 kpc) and his luminosity estimate for both com- 
ponents (2.4+J;9 kpc). 

The optical emission spectrum consists of a large num- 
ber of singly ionized metallic lines (Fen, Tin, Cm, Sill) 
and is dominated by the hydrogen lines of the Balmer series. 
V2028 Cyg possesses an extended gaseous and dusty enve- 
lope, which is indicated by several spectral feature s (IR ex- 



cess, forbidden lines, intensive Balmer emission). Marston 



& McCollum ( 2008 1 unsuccessfully tried to find evidence 
of larger-scale structures (lobes, nebula,. . . ) around the ob- 
ject using a narrow-band Ha filter. 



|Zickgraf k. Schulte-Ladbeck ( 1989[ ) assumed that a cir 



cumsteliar disc structure causes a direction-dependent po- 
larisation of the emitted photons. They determined an in- 
trinsic polarization of ~ 2% from their observations. 
On the basis of this value and the assumpti on that the scat- 
tering plane is parallel to the orbital plane, Zickgraf ( 2001 ) 
estimated an inclination angle > 30 — 45° . 

The gaseous and dusty envelope of V2028 Cyg makes 
it very difficult to determine the fundamental parame- 
ters of the object, since the commonly used stellar atmo- 
sphere models are not applicable in this case. The proper- 
ties of these objects can only be described by multidimen- 
sional models that combine time-dependent hydrodynamics 
and NLTE radiative transfer. This task is our ultimate goal. 



Our t ime-dependent hydrodynamical code (Votruba et al. 
20071, which is now extended into 2D, gives the tempera- 



ture and density distributions and v elocity of the material 
for the 2D radiative transfer code (Korcakova & Kubat 
2005 Korcakova et al. ]|2011 ). However, this is a long-term 
project complicated the role of different physical effects 
in the formation of the B[e] phenomenon remaining unclear. 
Therefore, we decided to focus our efforts here on the anal- 
ysis of the object's temporal behaviour that may reveal 
important physical mechanisms and narrow down the pos- 
sible system configurations. The first papers of the planned 
series focus on a description of the observed spectral prop- 
erties of selected B[c] stars over a time interval of several 
years. 

The observations, data reduction, and main spectral fea- 
tures are described in Sects. [2] and [3} The following Sect. [4] 
deals with the temporal variability of the emissions of Fe II 
(6318, 6456, 6417 A), [Oi] (6300 A), and Ha. In Sect. § we 
discuss the possible nature of V2028 Cyg based on our ob- 
servational results and give several criteria and suggestions 
for future modelling. 



2. Observations and data reduction 

The principal part of our data were obtained with the 2 m 
telescope located at the Ondfejov observatory, Czech 
Republic. We collected 88 spectra around the Ha line dur- 
ing 45 nights over nearly six years (December 1, 2004 - 
October 8, 2010). The wavelength interval in these spectra 
covers 6250 — 6770 A (Fig. 111. The resolving power is R ~ 
12500. Three additional spectral intervals were observed: 
7520 - 8030 A (seven spectra), 7700 - 8210 A (three spec- 
tra), and 8200 — 8700 A (three spectra). The Ondfejov spec- 
tra we re taken with the cl assical spectrograph at the Coude 
focus (ISkoda et al. ||2002h. 



This data set is supplemented by four spectra (ac- 
quired over the years 1994 — 1996) obtained with an echelle 
spectrograph ELODIE operated at the Cassegrain focus 
of the 1.93 m telescope at the Observatoire de Haute- 



Prove nce. The resolving power is 42000 (Moultaka et al. 
| |2004[ ) in the wavelength interval 4000 - 6800 A. 

Data reduction was performed using the IRAfTH data 



2004) 



reduction package, and the DCR program (Pych 
was used for cosmic-ray-hit elimination. 

3. Spectral features 

As mentioned above, the spectrum of V2028 Cyg contains 
both emission and absorption lines. The absorption spec- 
trum is dominated by the lines of the later component. 
The earlier component absorptions are observed in Hei 
lines (4026, 4144, 4472 A). The emission spectrum is domi- 
nated by Balmer lines (Ha - US) , which are about an order 
of magnitude stronger then the remaining emission lines 
in the spectrum. Most of the emission lines come from singly 
ionized metals (Fen, Tin, Cm, Sin). Forbidden lines are 
represented by [Fen] (4244, 4277, 4287, 7711 A) and [Oi] 
(6300, 636 4 A). The emission-line i dentification is mostly 
taken from iZickgraf & Stahl I (|1989h . 



3.1. Equivalent widths 

The equivalent widths of the lines were calculated by a nu- 
merical integration of the line profile using a trapezoi dal 



rule. The method of Vollmann & Eversberg (2006) is 



adopted for the error estimation. The signal-to-noise ratio 
(S/N) of individual spectra was calculated by performing 
a linear fitting of the continuum around the studied line. 
The weightin g function for the fit was taken from Mikulasck 



et al. 



([2003J). 

TEeS/N around the Ha is about 25 - 45. Higher S/N 
values cannot be achieved using the instrument that we 
used, as the maximum intensity of the Ha line is much 
stronger than that in the continuum. To obtain an un- 
saturated Ha line, the signal in the continuum cannot be 
high (300 — 900 ADU). The average spectrum was used 
when more then one observation was made during the night. 

The equivalent width measurements of the metallic lines 
are shown in Fig. [2] The mean values are -0.31±0.02 A(Fell 
6433 A), -0.90±0.05 A(Fen 6456 A), -0.35±0.03 A(Fen 
6417 A), and -0.49±0.04 A([0 1] 6300 A). The relative errors 



IRAF is distributed by the National Optical Astronomy 
Observatories, operated by the Association of Universities 
for Research in Astronomy, Inc., under contract to the National 
Science Foundation of the United States. 
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Fig. 1. Spectral interval around Ha observed at the Ondfejov observatory. Ha profile is cut to resolve other spectral 
lines. 
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Fig. 2. 
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of 10% — 26% do not, however, allow us to make a definite 
statement about the magnitude or periodicity of the varia- 
tions. 

In the case of the Ha line, the equivalent width vari- 
ations are clearly visible. The relative errors are around 
4%. Fig. [3] shows the changes in the Ha equivalent width 
and height. Both quantities exhibit very similar behaviours. 
Only data from the Ondfejov spectra are plotted in Fig. [3j 
therefore direct comparison is possible. 

3.2. Line profiles 

The metallic emission lines in the Ondfejov spectra 
are only slightly asymmetric, and their profiles are very 
close to a Gaussian. A higher resolution spectra from 
the ELODIE spectrograph clearly show the asymmetry 
of metallic emission lines (see Fig. |4| . 



Fig. 3. Equivalent widths (squares) and intensities (di- 
amonds) of the Ha line. Only Ondfejov data are pre- 
sented. The error bars of intensities are of the same size 
as the points. 



The blue side of the profile decreases more steeply then 
the red side, which is evidence of a mass outflow. This fea- 
ture is common in all metallic emission lines stronger than 
~ 1.8 times the continuum level. The line profile shape 
of the weaker lines cannot be precisely determined owing 
to the noise, but in some cases the same type of asymmetry 
is discernible, and also in the c ase of the [Oi] lines. Thi s 



feature was mentioned before by Zickgraf & Stahl ( 1989 



Balmer lines are strongly asymmetric and their asym- 
metry is somewhat the opposite of that of the metallic lines 
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Fig. 4. Emission line profiles of Fen 6456, 5317, 5018, 
and [O i] 6300. ELODIE spectrograph (night August 11/12, 
1995). 
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Fig. 5. Emission line profiles of Ha and H/3. ELODIE spec- 
trograph (night August 11/12, 1995). The line profiles are 
affected by the telluric lines in this figure. 



(Fig. [5]). The narrow peak tends to be shifted to the red side 
with respect to the lower parts of the core. The shift varies 
in time. The little depressions on the red side of the Ha 
represent the absorption of atmospheric water. 

The Ha line is the most pronounced line in the spectrum 
and its profile is strongly variable. The changes have a no- 
ticeable maximum and minimum and are seen not only 
in the line intensity and equivalent width (Fig. [3]), but also 
i n the b isector velocities. Previous observations of |Zickgrar| 
[] ( 2001 ) show even stronger variations then those present 
m our data set. He found that the intensity of the Ha line 
changed significantly over the years 1986 — 2000. 



3.2.1. Width and height of Ha 



The heights (the maximum flux /; Fig. [3]) and widths 
(Fig. [6]) of the Ha line vary during the observation run. 



The relative differences of the widths are as high as 24% 
and for the heights it is 21%. 

The absolute value of the equivalent widths \EW\ is 
correlated with height (correlation coefficient g — 0.75), 
but there is no correlation with the measurements of our 
so-called full widths (g = —0.06), which were measured 
at relative heights of from 0.1 (wings) to 0.9 (peak) with 
a step size of 0.05. To compare variations in the different 
relative heights, we divided each of them by its mean value. 
The relative variations in the individual heights were found 
to be very similar. 



1.2 
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Fig. 6. Variations in the full width of the Ha line. 
The mean value of the full widths measured in relative 
heights 0.1 — 0.9 is plotted in the graph. 



3.2.2. Ha bisectors 

A more detailed description of the profile variations was 
done using bisectors. Certain values of the bisectors were 
measured at the relative heights 0.1— 0.9 with a step of 0.05. 
We used the interpolation method of Steffen ( 1990| ) to ob- 
tain line profile points from the observed ones. 

Three parts of the Ha line can be resolved (Fig. [7]). 
The wings and the lower part of the core (0.1 — 0.4 of the rel- 
ative height) appear to have similar kinematical behaviours 
and in spite of the scatter, the curve with both its minimum 
and maximum is discernible in the plot. The minimum oc- 
curs on JD 2454130±40. The relative heights at 0.7 - 0.9 
represent the peak. The bisectors show that its radial ve- 
locity also slightly changes, but in a different way from 
the velocity variations of the wings. The minimum was 
observed on JD 2454540±30 here. The bisector velocity 
at half-maximum of the relative height of the line shows 
some scatter, but no noticeable trend or period. 



4. Radial velocities 

Since V2028 Cyg is assumed to be a binary, we investi- 
gated the radial velocity variations in its spectra. The ab- 
sorption lines of the cool component as well as the emission 
from the envelope were measured. 

The K type absorption lines are often blended and af- 
fected by noise in the Ondfcjov spectra. We therefore 
used the cross-correlation method for our measurements. 
The ELODIE spectrum reduced to the Ondfejov resolution 
serves as a template. Nine intervals without emission lines 
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Fig. 7. Bisector points of the Ha line in relative heights 
0.1 (wings) to 0.9 (peak). 



and atmospheric absorption lines were chosen for the mea- 
surements (6334-6344.5 A, 6398-6415 A, 6434-6440.5 A, 
6444 - 6453 A, 6607 - 6626 A, 6636 - 6664.5 A, 6680.5 - 
6693 A, 6694.5 - 6706.5 A, 6709 - 6720 A). The resulting 
radial velocities in Fig. [8] are mean values of these inter- 
vals and the errors are computed as the standard devia- 
tion in the arithmetic mean. The measured mean radial 
velocities vary slightly, but no periodical variations were 
found. To check the validity of the results we chose seven 
unblended absorption lines to measure radial velocities us- 
ing other methods (Fei 6355 A, Cai 6439 A, Crl 6630 A, 
Fe i 6647 A, Fe i 6648 A, Fe i 6704 A, Li i 6708 A) . The shifts 
in these lines were measured (Fig. [8]) using a) Gaussian 
fitting and b) profile mirroring. The profile mirroring was 
done using both automatic and manual procedures. The re- 
sults from both techniques agreed to within the error inter- 
vals. The lower panel of Fig. [8] shows the values obtained 
manually usi ng a program SPEFC[^] (developed by J. Horn, 
description in Skoda Il996). The scatter and error intervals 



are approximately double those of the cross-correlation re- 
sults, although the variations are very similar. There was 
a decline in the past 200 days and had been an indistinct 
maximum before it (Fig. . 

The radial velocities of metallic emission lines (Fig. [9]) 
were derived by Gaussian fitting. There is also some scatter, 
but this is primarily caused by the noise. No period can also 
be seen here. 

The radial velocities of the Ha line were measured 
for the peak as well as for the line wings. A polyno- 
mial fitting was used to measure the peak v elocit ies. We 
adopted the method of Mikulasek et al. | (|2003|) based 



We used a double-Gaussian method (Shafter 19831 
to measure the radial velocities of the wings. The method 
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Fig. 8. Radial velocities of the K type absorption lines 
measured by cross-correlation (upper panel) and profile 
mirroring (lower panel). 



Fell 6433 



[01] 6300 




3500 4000 4500 5000 5500 3500 4000 4500 5000 5500 



Fell 6456 



Fell 6417 




3500 4000 4500 5000 5500 3500 4000 4500 5000 5500 



Fig. 9. Radial velocities of metallic emission lines. 



consists of the determination of the line centre A by solving 
the equation 



S(A)K(X - A) dk = 0, 



(1) 



on the least squares method. The Ondfejov spectrograph where 
resolution provides typically six to eight points in the peak 
interval, which is sufficient to obtain a goo d fit. 



http://astro.troja.mff.cuni.cz/ftp/hec/SPEFO.263/ 



K(x) = cxp[-(x - a) 2 /2<r 2 } - exp[-(x + a) 2 /2a 2 } (2) 



and 5(A) is the investigated spectr um. The errors wer e 
computed according to the relation by Home et al. ( 1986 ). 
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Fig. 10. Radial velocities of the Ha peak (upper panel) 
and wings (lower panel). 



The parameter a adjusts the separation of the Gaussians 
and should be chosen according to the width of the wings 
and the core. Since the Gaussian parameter a is limited 
by the spectral resolution, the lowest possible value in our 
case is 0.25 A. However, a larger a value enables us to mea- 
sure the wider wavelength interval of the wings and it deals 
better with noise. We adopted values a = 5.0 A and a = 
1.0 A. The credibility of this method was verified by com- 
paring with the results obtained by the profile mirroring 
method. 

Fig. |10| (lower panel) shows the minimum and maximum 
of the raHIal velocities of the wings. The position of the min- 
imum is JD 2454155±17. The radial velocities of the Ha 
peak are plotted in the same figure (upper panel). The re- 
sults are different from those for the wing velocities. There 
is a large scatter in the data but it is possible to resolve 
the minimum at JD 2454683±33. 



5. Discussion and conclusions 

We have measured the equivalent widths and radial veloc- 
ities of several emission and absorption lines in the spec- 
tra of V2028 Cyg. We have also thoroughly investigated 
the profile variations in the Ha line. 

The emission lines of metals (the most common 
of which is Fen) have a characteristic bluewards skewed 
profile, which is indicative of an outflowing envelope. 
The line profile width at half-maximum is usually ~ 45 — 
55 km/s. The stronger lines (two times the continuum 
intensity) possess wings of widths up to ~ 200 km/s. 
Assuming spherical symmetry, this line profile is indica- 
tive of a slow (up to 100 km/s) outflow with a small ve- 
locity dispersion. Equivalent widths and radial velocities 
of the metallic emissions display some scatter but no peri- 



odic variations. The slight variations in the results may be 
ascribed to fluctuations in the envelope. 

The forbidden lines are in the chosen interval (A = 
6250 — 6770 A) represented mainly by [O i] emission. In ad- 
dition, their profiles are skewed bluewards but, however, 
lack the broader wings. They therefore probably originate 
in the outer, less dense layers of the envelope. 

The radial velocities of the K component measured 
for the single, most clearly defined absorption lines (de- 
rived using both the mirroring method and by fit- 
ting the Gaussian profile) have a similar behaviour 
to the Ha wings. However, the cross-correlation method, 
even if the chosen intervals are unaffected by water ab- 
sorption, gives less significant results. Since our spectra 
are affected by a large amount of noise, the radial veloc- 
ities obtained from the measurements of the most clearly 
defined spectral lines are probably more reliable than 
the latter. Furthermore, accurate observations are neces- 
sary to (dis)confirm the correlation of the radial veloci- 
ties of the Ha wings with those of the absorption lines. 
This relation can be important when determining the role 
of binarity in this system. The radial velocity variations 
of the absorption component in our data are relatively small 
(< 8 km /s) and fall into the same interval as the measure- 
ments of |Zickgraf (2001 1. This implies a long orbital period 
for V2028 Cyg with a lower limit of 25 years. 

The Balmer lines Ha, j3, and 7 are the strongest lines 
in the spectrum (approximately an order of magnitude 
stronger than the metallic emissions). All of them are 
asymmetric with the peak shifted redwards with respect 
to the wings. In most spectra, the lines bear a hump 
on the blue side, approximately at the height of the half- 
maximum. The H<5 emission is also clearly visible but any 
specific features of the profile are overlaid by the noise. 
From this set, only Ha has a sufficient coverage in the six- 
year-long data set from Ondfejov used to study its vari- 
ations. It is also a strongest emission line in the spec- 
trum. For this line, we measured the variations with time 
in the equivalent width, intensity, radial velocity, and bi- 
sector values. 

Our observations uncover some of the general proper- 
ties of the line profiles and correlations between measured 
quantities: 

- Ha line 

— The line profile is observed in nearly all spectra 
to be asymmetric. All of these observations show 
a red shift of the line peak relative to the centre de- 
fined by the wings. The symmetric line profile is very 
rare in our data. 

— The radial velocities and bisectors vary differently 
in the wings and peak. 

— The absolute value of the equivalent widthJ_EVF| is 
correlated with the maximum flux I (Fig. uty . 

— The absolute value of \EW\ is anti-correlated 
with the radial velocity of the wings v rw (Fig. [3] 
and [10]). 

- The equivalent width of [Oi] 6300 A, Fen 6427, 6433, 
and 6456 A lines do not show any significant changes. 

- The radial velocities of the K absorption component 
may correlate with the Ha wings. However, further ob- 
servations of higher resolution are necessary to confirm 
this property. 
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These results can be compared with the results of future 
modelling. In particular, the anticorrelation of the absolute 
value of the Ha equivalent width and the radial velocities 
of its wings probably strongly restricts on the proposed 
model. 

V2028 Cyg is a very complicated system to model. It is 
necessary to take into account the binary nature of the ob- 
ject. A particularly large problem is the extended circum- 
stellar region, which probably forms a disc-like structure. 
These properties require a combination of at least 2D 
time-dependent hydrodynamic and NLTE radiative trans- 
fer in moving media. Although our observations have not re- 
vealed the real nature of V2028 Cyg, they specify the model 
properties for future investigation. Particular model prop- 
erties can be tested using our observations. 

Acknowledgements. We thank especially Anatoly Miroshnichenko 
for his valuable and inspiring advice and remarks. We also thank 
P. Hadrava, P. Koubsky, and P. Harmancc for their advice and Adcla 
Kawka for taking some of the spectra in her observation time. This 
research is supported by grants 205/09/P476 (GA CR), 205/08/H005 
(GA CR), MUNI/A/0968/2009, and MSM0021620860 (MSMT CR). 
The Astronomical Institute Ondfcjov is supported by a project 
AVOZ 10030501. 



References 

Allen, D. A. 1973, MSRSL, 5, 335 

Allen, D. A. 1974, MNRAS, 168, 1 

Allen, D. A., & Swings, J. P. 1976, A&A, 47, 293 

Arkhipova, V. P., & Ipatov, A. P. 1982, SvAL, 8, 298 

Bergner, Yu. K., Miroshnichenko, A. S., Yudin, R. V., ct al. 1995, 

A&AS, 112, 221 
Castor, J. I., Abbott, D. C, & Klein, R. I. 1975, ApJ, 195, 157 
Corporon, P., & Lagrange, A.-M. 1999, A&ASS, 136, 429 
Home, K., Wade, R. A., & Szkody, P. 1986, MNRAS, 219, 791 
Korcakova, D., & Kubat, J. 2005, A&A, 440, 715 
Korcakova, D., Nagel, T., Werner, K., Suleimanov, V., & Votruba, V. 

2011, A&A 529A, 119 
Lamers, H. J. G. L. M., Zickgraf, F.-J., de Winter, D., Houziaux, L., 

& Zorec, J. 1998, A&A, 340, 117 
Marston, A. P., & McCollum, B. 2008, A&A, 477, 193 
Merrill, P. W., Burwell, C. G., & Miller, W. C. 1942, ApJ, 96, 15 
Mikulasek, Z., Ziznovsky, J., Zverko, J., & Polosukhina, N. S. 2003, 

Contrib. Astron. Obs. Skalnate Pleso, 33, 29 
Miroshnichenko, A. S. 2007, ApJ, 667, 497 

Miroshnichenko, A. S., Chentsov, E. L., & Klochkova, V. G., et al. 

2000, A&AS, 147, 5 
Moultaka, L., Ilovaisky, S. A., Prugniel, P., & Soubiran, C. 2004, 

PASP, 116, 693 
Pych, W. 2004, PASP, 116, 148 
Shatter, A. W. 1983, ApJ, 267, 222 
Skoda, P. 1996, ASPC, 101, 187 

Skoda, P., Slechta, M., & Honsa, J. 2002, AICz 90, 3 

Steffen, M. 1990, A&A, 239, 443 

Vollmann, K., & Eversberg, T. 2006, AN, 9, 865 

Votruba, V., Feldmeier, A., Kubat, J., & Ratzcl, D. 2007, A&A, 474, 
549 

Zickgraf, F.-J., 2001, A&A, 375, 122 

Zickgraf, F.-J., & Schulte-Ladbeck, R. E. 1989, A&A, 214, 274 
Zickgraf, F.-J., & Stahl, O. 1989, A&A, 223, 165 

Zickgraf, F.-J., Wolf, B., Stahl, O., Leitherer, C, & Klare, G. 1985, 
A&A, 143, 421 



